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ABSTRACT 

A study was conducted to determine what the parameter s of an acoustic liner should 
be to maximize liner damping. This was accomplished by solving the wave equation ex- 
actly for the frequency and damping of low order tangential and radial modes in a lined 
cylindrical cavity. The results, presented as dimensionless frequency and damping for 
each mode plotted in the impedance plane, indicate that a, the acoustic absorption coef- 
ficient, can be used to design a liner provided that the liner real impedance is larger 
than unity. 
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ON THE DESIGN OF ACOUSTIC LINERS FOR ROCKET ENGINES: 

MAXIMUM DAMPING AS A DESIGN OBJECTIVE 
by Bert Phillips 
Lewis Research Center 

SUMMARY 

A study was conducted to determine what the parameters of an acoustic liner should 
be to maximize liner damping. This was accomplished by solving the wave equation ex- 
actly for the frequency and damping of low order tangential and radial modes in a lined 
cylindrical cavity. The results, presented as dimensionless frequency and damping for 
each mode plotted in the impedance plane, indicate that a, the acoustic absorption coef- 
ficient, can be used to design a liner provided that the liner real impedance is larger than 
unity. 


INTRODUCTION 


Extensive testing at the Lewis Research Center and elsewhere has shown that an ar- 
ray of Helmholtz resonators can, under certain conditions, damp combustion instability 
in rocket engines (refs. 1 and 2). The conditional success of some liners is indicated by 
the lack of correlation between high theoretical absorption coefficient, a, and liner ef- 
fectiveness as shown in reference 1, page 30, and reference 3. A possible reason for the 
disparity between theory and experiment is that, according to reference 4, the damping of 
acoustic waves in a cavity (combustion chamber) is maximized by maximizing a only 
under the condition that the wave frequencies are considerably higher than the fundamen- 
tal cavity or chamber resonance frequencies. The results of reference 5, however, in- 
dicate that the liner, if it is to be effective, must absorb frequencies near to the lowest 
resonant frequencies of the chamber. Consequently, it is not clear that maximizing the 
acoustic absorption coefficient, a , is the proper design objective. 

To determine what the parameters of an acoustic liner should be to maximize the 
damping of a liner, the exact solutions for the tangential and radial acoustic modes in a 
closed end (no flow) cylindrical cavity with absorbing walls were obtained. The approach 



is similar to that presented in reference 6 without the combustion driving and nozzle 
damping. The study presented in reference 6, however, made no attempt to investigate 
what parameter of a liner would maximize the liner damping. To facilitate the use of the 
results by engine designers, the results are presented as the normalized frequency and 
damping for each mode as functions of the liner parameter 0, the real part of the acoustic 
impedance (resistance) and X, the imaginary part of the acoustic impedance (reactance). 
The solutions were obtained by the use of a program for calculating complex Bessel func- 
tions reported in reference 7. Since combustion instability only involves the lower order 
resonant modes, only the lower order solutions will be presented. 


THEORY 


The equation for acoustic waves in a cylindrical cavity with closed ends, no mean 
flow, and uniform gas properties is: 
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where 
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p oscillatory pressure, lbf/ft 

r radial distance, ft 

z axial distance, ft 

(p azimuthal angle, rad 

c cavity sonic velocity, ft/sec 

t time, sec 

For a transverse mode with no axial dependence and exponential time variation, a solu- 
tion for equation (1) is: 
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J Bessel function of first kind of order n, n = 0, 1 . . . 

n 7 

a Q cavity radius, ft 

o> angular frequency, rad/sec 

2 



k temporal damping constant, 1/sec 

P n ,k n distribution coefficients resulting from cavity walls that are not infinitely hard 

9 

P maximum value of oscillatory pressure, Ibf/ft 

Substituting equation (2) into equation (1): 
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This expression gives the relation between the oscillation frequency, the decay con- 
stant, and the distribution coefficients. The next step is to solve for the distribution co- 
efficients in terms of the wall resistance and reactance, 9 and X. 

The radial acoustic particle velocity into the wall is: 
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u r radial velocity, ft/sec 

3 

p gas density, lbm/ft 

2 

g gravitational constant, (lbm/lbf)(ft/sec ) 

c 

Substituting equation (2) into equation (4) and setting n = 1 (first tangential mode) 
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For clarity, the quantity (p n + ik n ) is redefined, as in reference 7, as Re , where R 
is the distribution coefficient magnitude and h the distribution coefficient angle. 


7 r(p n + ik n ) = R cos(h) + i R sin(h) 


( 6 ) 


Dividing equation (5) into equation (2) and defining r] = coa Q /7rc the following equation re- 
sults: 
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The quantity P/[(p/g c )u r c] is the impedance which can also be expressed as the sum of a 
real and imaginary part, 
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where 

rj wavelength ratio defined by o>a o /7rc = D/X 

D cavity diameter 

X wavelength 

0, X normalized resistance and reactance 

By choosing values of R and h, equations (7) and (8) can be solved for 0/rj and X/17. 
Combining equations (3) and (6), it can be shown that 

kci 

— - = R sin(h) (9) 
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These are the normalized dampii^ rates and resonant oscillation frequencies of the cavity 
for the given values of R and h. Therefore, by specifying R and h, wa Q /c, and r) 
can be calculated and, from equation (8), 0 and X can be determined. 


RESULTS AND DISCUSSION 
Defining Modes of Oscillation 

Modes of acoustic wave oscillation in cavities are generally defined by the values of 


4 



wa Q /c that satisfy the condition that u r , the radial velocity at the wall, is zero. For the 
case of nonzero admittance, and, in particular, for 9 = 1, mode definition becomes quite 
difficult. The approach taken in the present report was to define the modes based on the 
value of R, defined by equation (6). This approach has some justification since, in the 
limit of hard walls (h = 0°, 9 = «), R = coa Q /c and each mode has a single value of R. 

To determine the behavior of R in the region of interest ("soft walls’ 1 ), lines of constant 
R were plotted in the 9 - X plane, as shown in figure 1. The Jq Bessel function was 
chosen for illustration. For all of the R values shown in the figure, R at h = 0° is on 
the abscissa (9 = 0 ); as h is increased, the line spirals in toward the origin until, at 
h = 90°, all values of R intersect the origin. For R = 0. 5, the h = 0° point is at 
X= 3.074, 0 = 0°. As R is increased, the h = 0° point on the abscissa moves left to- 
ward the origin until, at R= 2.42, the line intersects itself at h = 0° and 90°. Between 
R = 2. 42 and 3. 05, the lines of constant R cross, An example is R = 2. 80 which 
crosses at point A on figure 1, corresponding to h = 3° and 66°. Continuing to increase 
R beyond 3. 05 moves the h = 0° point to the left until, at R very close to 3. 83, the 
h = 0° point is at X = -°°. At R exactly equal to 3. 83, the h = 0° point is at 9 = °°, 

X = 0. For R > 3. 83, the h = 0° point swings to X = +°° and again moves to the left 
with increasing values at R. The process continues as R increases, with lines of con- 
stant R crossing within certain ranges. 

The same results are obtained with Jg, etc. , giving rise to a series of curves in 



Figure 1. - Lines of constant distribution coefficient magnitude as function of 
real and imaginary wall impedance. 



Figure 2. - Values of distribution coefficient magnitude 
and angle designating mode boundaries. 


the R, h plane defined by the values of R and h at the crossings. Four of these curves 
are shown in figure 2, On the Jqq - Jq^ curve, the two points, labeled A, correspond 
to point A on figure 1. These curves in the R - h plane are arbitrarily taken as the 
boundaries of modes. All values of R and h below the Jqq - Jqj curve are in the Jqq 
mode, all values of R and h between the Jjo “ ^11 an( *the **n " ^12 curves are in 
the J|j mode, etc. Thus by defining the value of R and h, we can confine out discus- 
sion to a particular mode. 


Modal Resonant Frequencies and Damping 

Using the previously obtained R-h criterion, the fundamental dimensionless reso- 
nance frequencies and damping of each mode can be calculated as functions of 9 and X. 
These results are displayed as lines of constant a>a Q /c and ka Q /c in the 0,X plane. 
The results for the Jqq mode are presented in figures 3(a) and (b). The Jqq mode is 
presented first although it has no real physical significance since the fundamental reso- 
nance frequency of a Jqq mode for hard walls is 0. 

The lines of constant a)a Q /c in figure 3(a) resemble those of constant R. In the re- 
gion < 0 ^ ® >, all the lines of constant frequency converge, resulting in a zone 

[p >X > -0. 6J 
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Figure 3. - Lines of constant frequency and constant damping as function of real and imaginary wall impedance for Jqq mode. 


where the resonant frequency is not well defined, that is, small changes in either 9 or 
X can result in large changes in coa Q /c. As Wa 0 /c is increased to 2. 42 and greater, 
the lines of constant frequency intersect in the same manner as lines of constant R in 
figure 1. The interpretation can be made that, in this region, both the Jqq and the Jqj 
mode have the same resonant frequencies. The results indicate that the resonant fre- 
quency for the mode can be specified readily anywhere in the 6 - X plane except for the 
region already mentioned ( 9 < 0. 8). 

The behavior of the damping for the Jqq mode was determined in the same manner 
as the frequency and is shown in figure 3(b). The lines of constant damping begin at the 
origin and spiral clockwise until the corresponding values of R and h on the lines of 
constant ka Q /c correspond to the upper bound of the mode, where the line is terminated. 
The region of maximum damping, indicated by the cross-hatching, on the figure corre- 
sponds to the upper bound of the mode, and, in addition, corresponds to the region of ill- 
defined frequency of figure 3(a). Thus, both frequency and damping for the mode can be 
specified except for the same cross hatched regions. Although the discussion of the 
damping characteristics of the J Q0 mode is, by no means, completed, it may be more 
instructive to go on to the next solution of Jq, the Jqj mode. 

The lines of constant frequency and damping are presented on figures 4(a) and (b). 
The results on figure 4(a) are similar to those presented in figure 3(a), however, the 
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Figure 4. - Lines of constant frequency and constant damping as function of real and imaginary wall impedance for J Q1 mode. 


resonant frequency of the mode in the limit of hard walls is not zero. The line rep- 
resenting coa Q /c = 3. 83 continues upwards, representing the resonant frequency of the 
J Q1 or first radial mode for hard walls. In the limit of large 0, as the reactance, X, 
becomes negative, the line of constant frequency for values less than 3. 83 fall to the left, 
that is, the resonant frequency decreases. If there is a positive wall impedance, X > 0, 
the Jqj resonant frequency becomes greater than 3. 83. This result confirms the ex- 
perimental and theoretical approximation results of reference 9. In the same manner as 
figure 3(a), there is a region of poor definition in the frequency where the lines of con- 
stant ct>a 0 /c converge. 

The lines of constant damping for the Jq^ mode are presented in figure 4(b). The 
results are quite different from those of figure 3(b) and show lines of damping that are ap- 
proximately closed concentric curves except for the cross-hatched area. If the high 
damping in the cross-hatched area can, for the moment, be disregarded, the results may 
be interpreted as defining zones of maximum damping in the region of 0-1, X = 0. For 
comparison, lines of constant a , the absorption coefficient, are plotted in the impedance 
plane and are shown in figure 5. The absorption coefficient, a, which maximizes in the 
same region, 9 = 1, X = 0, is an appropriate design criterion, since maximizing a 
would maximize the true damping. 
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Figure 5, - Lines of constant absorption 
coefficient as function of real and 
imaginary wall impedance. 


It is important, however, to investigate any advantage in designing a liner for the re- 
gion of maximum damping indicated by the cross-hatched area. Consider the Jq ^ mode 
for 9 less than unity with X = -2; if the liner tuning could be gradually improved so as 
to make X less negative, the resonant frequency would gradually decrease as per fig- 
ure 4(a) until the common line representing all the frequencies was reached. 

As X approaches -0. 4 from the left, the damping, according to 4(b), sharply in- 
creases at the cross-hatched region. The frequency, however, as per figure 4(a), is 
near the common line and the system is now operating in the regime where both the Jqq 
and Jq j modes have the same frequency. Any system, given a situation where it can 
resonate at the same frequency with two widely different damping rates will appear to de- 
cay at the lower damping rate since the higher damping wave will not be evident. In the 
case under discussion, the system will decay at the lower damping rate corresponding to 
the Jqq mode. In addition, the oscillatory pressure profile will generally shift to that 
corresponding to the mode that has the lower damping. Consequently, it does not appear 
that any advantage could be gained by designing for the cross-hatched region of damping. 

If, therefore, the cross-hatched region were disregarded in designing a liner, a could 
properly be used for the first radial mode, provided that the values of 9 were maintained 
above unity, and away from the cross-hatched region. 

The next system to be investigated was the J 10 or first tangential mode. A plot of 
the lines of constant frequency in the 9 - X plane is shown in figure 6(a). The results 
appear to be similar to those presented in figure 4(a). In this case, however, the reso- 
nant frequency corresponding to the first tangential mode with ’’hard walls” (0 ^ °°), is 
1. 84. In addition, the effects of positive and negative impedance on the resonant fre- 
quency in the region of large 9 are consistent with those of figure 4(a). The region cor- 
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responding to ill-defined resonant frequency is an area in the same position as that region 
for the Jq modes. 

The lines of constant damping for the J^q mode are shown on figure 6(b). The re- 
sults resemble those for the J Q q mode as shown in figure 3(b) with the region of maxi- 
mum damping in the same position. Following the same approach as with the Jq modes, 
the discussion of the J.q mode is now in order. The frequency and damping results for 
the Jjj mode are presented in figures 7(a) and (b). The frequency in the limit of 9 - °° 
for this mode is 5. 33 corresponding to the first tangential-first radial mode with ’’hard 
walls. " The damping for the mode is, with the exception of the cross-hatched area, 
a set of closed curves of lines of constant damping with the maximum values in the 9 ~ 1, 
X ~ 0 area. 

Continuing both the Jq and J^ solutions, a set of figures for both the frequency and 
damping can be obtained. In all of the cases, the lines of constant damping will be closed 
curves with the exception of the cross-hatched region which is common to all the modes. 
This cross-hatched region can be regarded as the operating point at which modes can 
readily be altered at constant frequency. Given a system with a driver with finite band- 
width, operation anywhere in the region 9 < 1, X < 0 will encourage the system to 
switch to the lower damping mode. 


CONCLUDING REMARKS 

The conclusion that one can come to with regard to the results presented is that a 
can be used to design a liner which will damp low order transverse cylindrical modes 
provided that the real impedance is maintained at values greater than unity. Operating at 
a value less than unity with negative imaginary impedance may result in the damping of 
one mode and the appearance of another. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, October 9, 1968, 

128-31-06-05-22. 
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